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Abstract. - The electron tunneling through a molecular junction modeled by a single site weakly 
coupled to two leads is studied in the presence of a time-dependent external field using a master 
equation approach. In the case of small bias voltages and high carrier frequencies of the external 
field, we observe the phenomenon of coherent destruction of tunneling, i.e. the current through the 
molecular junction vanishes completely for certain parameters of the external field. In previous 
studies the tunneling within isolated and open multi-site systems was suppressed; it is shown 
here that the tunneling between a single site and electronic reservoirs, i.e. the leads, can be 
suppressed as well. For larger bias voltages the current does not vanish any more since further 
tunneling channels participate in the electron conduction and we also observe photon-assisted 
tunneling which leads to steps in the current-voltage characteristics. The described phenomena 
are demonstrated not only for monochromatic fields but also for laser pulses and therefore could 
be used for ultrafast optical switching of the current through molecular junctions. 
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Introduction. — Electronic transport through molec- 
ular wires and junctions has recently attracted much at- 
tention experimentally as well as theoretically [1-3]. Un- 
der the influence of a bias voltage and because of the 
coupling to the leads which act as electron source and 
drain, a current through the molecular junction is estab- 
lished. When an external time-dependent field, such as 
d laser field or an additional ac voltage is applied to the 
molecular junction, several interesting effects arise. One 
jihenomenon is the well-known photon-assisted tunneling 
(PAT) [4] . It was studied already in the early 1960's exper- 
imentally by Dayem and Martin [5] and theoretically by 
Tien and Gordon using a simple theory which captures al- 
ready the main physics of PAT [6] . The main idea is that 
an external field periodic in time with frequency to can 
induce inelastic tunneling events when the electrons ex- 
change energy quanta fku with the external field. Another 
important effect is the famous phenomenon of coherent 
destruction of tunneling (CDT). Grossmann et al. [7-9] 
first studied this effect and found that tunneling can be 
quenched in a periodically driven quantum system. In 
the context of molecular wires, this phenomenon can be 
explained using Floquet theory in the case of a periodic 
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laser field [10], and CDT occurs for certain amplitudes 
of the laser field at fixed frequencies [11-15]. Different 
scenarios of controlling the tunneling in molecular wires 
and quantum dots have been proposed based on different 
mechanisms [4,14-18]. Also current-induced light emis- 
sion in molecular junctions has been studied [19]. 

In the current paper we focus on the tunneling through 
a single-site molecular junction. This might be a quan- 
tum dot (though the temperatures in the current examples 
are rather high for quantum dots) or a single molecular 
level acting as a molecular wire. The theoretical founda- 
tion is a density matrix formalism using a perturbative 
treatment within the molecule-lead coupling to second or- 
der. Applying this technique it is possible to calculate the 
time-dependent population in and the current through the 
molecular junction under the influence of a time-varying 
external fleld [15,20]. Since the effect of the external fleld 
on the coupling between molecule and leads is treated ex- 
actly and not neglected as for example in Redfleld the- 
ory, effects based on the influence of the laser on this 
coupling can be investigated. The tunneling between the 
molecule and the leads can, for small bias voltages, be sup- 
pressed by a monochromatic laser. This result is different 
from previous studies [13-15,21] in which the current van- 
ished because of the CDT between the sites of the wire. 
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In Ref. [13] the possibility of CDT for a single site was 
briefly mentioned as a limiting case of a wire with two 
sites having a large intersite coupling, but it was not fur- 
ther explored. Additionally we demonstrate that not only 
monochromatic laser fields but also laser pulses can lead 
to CDT. This opens opportunities for building optical cur- 
rent switches with a time resolution on the femtosecond 
scale. h=l is used throughout the rest of this Letter. 

Model. — The total system Hamiltonian H{t) = 
Hs{t) + Hb + HsR includes three parts: the relevant 
system Hs{t) mimicking the molecule, the reservoirs 
modeling the two leads and the system-reservoir coupling 
HsR. The creation and annihilation operators of electrons 
with spin a are denoted by and Ca, respectively, so that 
the description of a single-site molecule reads 



The excited state in the investigations by Bruder [23,24], 
Oosterkamp [25-27] and others has the same effect [28,29]. 

Since in most cases one is only interested in the time evo- 
lution of the relevant system, i.e. in this case the molecule, 
a time-local quantum master equation (QME) based on 
a second-order perturbation theory in the molecule-lead 
coupling has been developed for the reduced density ma- 
trix of the molecule ps{t) [15,20] 



dpsjt) 
dt 



-iCs{t)ps{t) 



(3) 



^[K^„, Axx'a{t)ps{t) - ps{t)Axx'a{t)] 



with auxiliary operators for the molecule-lead coupling 

Hs{t) = y2{so-fiE{t))cic, + Uc\c^c\ci (1) A,,,„{t) = f dt'C,{t - t')Us{t,t')K,,,, (4) 

" t 

A^Mt) = [ dt'C*,,{t - t')Us{t, t')K^., . (5) 
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with spin-independent on-site energy eg and electron in- 
teraction U within the doubly occupied states. The time- 
dependent term —iiE{t) describes the effect of the external 
field E{t) and for simplicity we assume that the propor- 
tionality factor —p equals unity. The external field is of 
the form E{t) = Eo{t) coswt with a possibly time- varying 
amplitude Eq (t) . The two leads coupled to the molecular 
junction are mimicked as electron reservoirs in thermal 
equihbrium by Hr = Y.q^^qab\„hqa- Here 6^^ and bqa 
denote the creation and annihilation operators of an elec- 
tron with spin a in reservoir mode LOq. Due to the assumed 
thermal equilibrium of the electronic leads, the occupation 
expectation values of the reservoir modes are determined 
by (b^q^bq'cr) = upiijJq — EF)Sqq' where Up is the Fermi 
function and Ep the Fermi energy. In further derivations 
we will only refer to the left lead but the formalism has 
to be applied to the right lead coupled to the wire as well. 
The coupling of the left lead to the molecule is given by 



Here Cs{t) = [^^(t),*] is the Liouville operator, 
Us{t,t') = T+exp{— i J^, drCsiT)} the time evolution op- 
erator and Cxx'{t) = trii{e'-'^«*$xe~'^^*$a;'PR} the reser- 
voir correlation functions with spin-independent reservoir 
density matrix pR. T_|_ is the time-ordering operator. 

The properties of the Fermionic reservoirs are described 
by a single quantity, the spectral density Jr{lu). For 
a dense spectrum, Jr{lo) is a smooth function and one 
can approximate it by a numerical decomposition into few 
Lorentzian functions [15] 
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K.a^xa = Y.^Vqc\^bq, + V*blc^,) (2) 
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with = = Y.qyqbqa, Ki„ = K^^ = , and a 
wire-lead coupling strength Vq. For the coupling of the 
molecule to the right lead similar equations hold. 

The calculations are performed with and without elec- 
tron interaction in the molecule. Neglecting electron in- 
teraction all electron-carrying states of the molecule are 
degenerate with energy Sq. In the case of electron inter- 
action U, double occupancy leads to a state with energy 
2£o + U. When the bias voltage is small and sq is between 
the Fermi energies Ep ,. and Epj of the right and the left 
lead, respectively, only the tunneling through the channel 
with energy eq leads a current, Iq. When the value of 
the bias voltage is above U, an extra channel opens and 
there is a shoulder in the current. Thielmann et al. [22] 
used these energy levels to explain the shoulders in the 
tunneling current for a quantum dot at low temperature. 



with real fitting parameters pk, and F^. Using the 
theorem of residues and denoting the Fermi function as 
np yields 

m 



fe=i 



2i 



m-\-m' 



^ fc=l 
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m-\-m' 



p721* 



(8) 
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with the abbreviation = Ofc ± iVk and the Matsubara 
frequencies Vk = i ^Trfc+Tr ^ rjnj^^ infinite sums over the 
Vk can be truncated at a finite value depending on the tem- 
perature T and the spectral width of Jr{u). With these 
forms of Ci2 and C21 one can obtain a set of differential 
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equations for the auxiliary density operators 

+iL'^L'Jt), (9) 
+ hLy^L'At) (10) 

with A^^,At) = ET=r' K.'At) and A^^,At) = 
^m+m As has been detailed previously [15], the 

coupled differential Eqs. (3, 9, 10) can now be numerically 
integrated applying, e.g., the Runge-Kutta scheme. 

Using the electron number operator of the left lead with 
the summation performed over the reservoir degrees of 
freedom Ni = '}2qcr^\Aq<^^ expression for the current 
is given by [15] 

hit) = e^^iv{Nips{t)} = -ieir{[NuH(t)]ps[t)} 

= 2eRe (trg {c\^Ai2a{t)ps{t) - c\^ps{t)li2a{t)]) (H) 

Here e denotes the elementary charge. This equation 
describes the current Ii{t) from the left lead into the 
molecule. A similar expression holds for Ir{t) from the 
right lead into the molecule. In a steady state and after av- 
eraging over one period of the driving field, // and 1^ have 
the same magnitude but opposite signs and therefore a to- 
tal transient current through the molecular junction can 
be defined as I{t) = {Ii{t) — Ir{t))/2. The time-dependent 
average current I will be determined below by averaging 
I{t) over five periods of the highly oscillating carrier field. 

A simple spectral density with only one Lorentzian 
{m—1) was chosen. With ili = eq we locate the maximum 
of the coupling at the site energy. By using Fi = buj, the 
coupling between the leads and the system is almost in 
the wide-band limit. Choosing pi = OMujUiT^ we obtain 
a maximum coupling strength of 0.01 lu which is smaller 
than the thermal energy fe^T — 0.025aj and much smaller 
than the external field energy uj. So the results below are 
all within the high-frequency limit. 

According to the approximate Tien-Gordon model [4, 
6] for monochromatic external fields which is based on 
a simplified scattering picture, the rectified dc currents 
through ac-driven molecular junctions are determined as 
[4] 

Itg= E JU^)llieo + n.) = J2l„ (12) 

n— — oo ^ ^ n— 

where the current in the driven system is expressed by 
a sum over contributions of the current /°c(^o + nu) in 
the un-driven case but evaluated at side-band energies 
£o + nuj shifted by integer multiples of the photon quan- 
tum and weighted with squares of Bessel functions. Note 
that the partial currents /„ contain contributions from 
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Fig. 1: Current induced by a monochromatic laser with am- 
plitude i?o = 2.405a). Average current /, Tien-Gordon current 
Itg, and partial currents /o to /s vs Vb/2aj for fcsr=0.025ci;, 
[7=0. From top to bottom: Itg, 1, li, I2, I3, lo- Shoulders in 
the currents are at Vb/(2uj)=l, 2, 3, etc. 

±n. The term J„(i?o/w) denotes the nth-order Bessel 
function of the first kind. The photon absorption {n > 0) 
and emission (n < 0) processes can be viewed as creat- 
ing effective electron densities at energies Sq ± nuj with 
probability J'^{Eq/uj). Already Dayem and Martin [5] 
showed clear steps in the tunneling current induced by 
these two processes. Although "a great deal of information 
can be extracted from simple models like the Tien-Gordon 
model" [4] in many cases one needs more sophisticated ap- 
proaches. In double quantum dot theories, for example, 
the Tien-Gordon formula has been recovered in the limit 
of weak inter-dot coupling [4, 30-32] . 

In the following we use a bias symmetric with respect 
to £0, i.e. Ef,i = £0 + Vbl2 and Ep.r = eo - Vbl2. The 
current Itg obtained with the help of the Tien-Gordon 
theory is compared to the average current I{t) obtained 
from the QME described above. The partial currents /„ 
defined in Eq. (12) are used to qualitatively explain the 
steps in the I-V characteristics. In the results shown below 
we first restrict ourselves to a monochromatic laser field 
with constant field amplitude of 2.405 lo, i.e. a zero of the 
zeroth-order Bessel function Jci{Eo/uj = 2.405)=0. 

Monochromatic laser, U — 0. In this case, en- 
ergy levels are equal to £0 regardless whether there is al- 
ready an electron with opposite spin on the site or not, 
and lead to the partial current Iq . The PAT-induced states 
have also equal energies £0 ± nuj leading to the inelastic 
current contributions /„. In Fig. 1 the average current I 
shows a step at 14/2 = oj. When Vb/2 is smaller than w, 
the current vanishes. At these low bias voltages channels 
with n ^ cannot participate in the conduction since the 
energies of the PAT-induced states are not in the window 
between the Fermi energies of the right and left contacts. 
Because Eq/uj — 2.405, Iq also vanishes due to CDT. For 
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other ratios of Eq/lu it is non-zero. When Vh/2 — oj, the 
left Fermi energy Epj — Eq + oj equals the PAT-induced 
energy and Ep^r = £o~^- Therefore the current Ii jumps 
to a finite value, so do Itg and /. A comparison of Itg 
and / in Fig. 1 shows that each step in these I-V curves is 
related to an inelastic current /„ . The slope of the steps is 
of course directly connected to the temperature. The sum 
Itg of all partial currents for n=0,. . ., 3 agrees reasonably 
with the current / from the QME, especially for small bias 
voltages. 

For the current through two coupled quantum dots the 
Tien-Gordon result is perturbative in the interdot cou- 
pling, as mentioned above [30-32]. For that case, as in 
our study, the Tien-Gordon approximation overestimates 
the current /„. The deviations between the average cur- 
rent / and the Tien-Gordon current Itg become larger 
with increasing bias. In the present study and for a small 
bias voltage the average current stems only from the main 
contribution Iq and for this case the Tien-Gordon results 
match the present results for different laser amplitudes 
Eq (not shown). For larger bias voltages, the energetic 
positions of the steps are equal though the Tien-Gordon 
approach seems to overestimate the contributions from the 
PAT-induced states sq ±nuj. 

At first sight it might be astonishing that CDT also 
works for one site coupled to reservoirs since the infinitely 
large reservoirs could be expected to destroy any coherent 
effect. Let us go back to the two- level system [7, 9] in 
which it does not matter whether both or just one level is 
periodically driven to be able to observe CDT. It is also 
not important if the static level is replaced by a continuum 
of states, since the driving only modifies the coupling. As 
can be seen from the present results the situation is even 
not changed if the continuum of states is replaced by a 
reservoir or if a second one is added. In all these situations 
the coupling between the sites can be suppressed if the 



correct amplitude of the external field is applied. 

Monochromatic laser, < U < u). Let us first 
concentrate on this case of small Coulomb interaction. 
When the bias voltage Vb is small, the current is estab- 
lished through the channel with energy Eq- As soon as 
Vb/2 approaches uj, an extra channel with energy Eq + U 
is opened by the electron interaction and a step in the 
current-voltage characteristics appears. Under the influ- 
ence of the external field, the PAT induces additional 
current channels, i.e. states with energies Eq ± nu and 
Eq + U ± nuj. As mentioned above, the bias is increased 
symmetrically with respect to eo, i.e. Ep.i = eq + Vb/2 
and Ep^r — eq — Vb/2. So a step in the I-Vb curve appears 
whenever an additional conduction channel can partici- 
pate in the tunneling. Therefore one finds steps at Vb/2 
values nuj — J7, nuj and nu -\- U. In Fig. 2 with U =0.4 lu 
one observes respective shoulders at Vb/{2cj)—0.6, 1.0, 1.4, 
1.6, 2.0, etc. The CDT phenomenon occurs again for small 
bias voltages and disappears at Vb/2 = 0.6a;. The partial 
current Iq has, in principle, a shoulder at Vb/2 = OAuj but 
it is not visible in Fig. 2 because Iq = due to CDT. 

Monochromatic laser, U > u). In this case the 
current-voltage curve behaves similar to that for < U < 
uj. Current steps still appear at Vb/2 values nu and noj + U 
when the energy levels match the higher (left) Fermi level. 
The energy levels Eq ~\- U ~ nuj of some conduction chan- 
nels can be larger than Eq and might coincide with the 
higher Fermi level, but the states with Eq — nut have again 
to be compared to the lower (right) Fermi energy. As 
shown in Fig. 3 for [7=1.4 w one observes the first shoul- 
der at Vf,/(2a;)=0.4 and this step relates to the channel 
energy Eq + U — uj. Further shoulders can be observed at 
T4/(2w)=0.6, 1.0, 1.4, 1.6, 2.0, etc. All energetic step po- 
sitions except the one at Vb/(2u;)=0.4 are the same as in 
Fig. 2 due to the fact that the electron interaction strength 
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Fig. 4: Current induced by laser with Gaussian shape of the 
envelope with amplitude Eq — 2.405tJ exp (— — fo)^/(2(j^)) 
with to=400 fs and a=80 is. Average current I vs time for 
kBT=0.025w. The two top lines are for a bias 14/2=2.4tj. The 
two lower lines are for a small bias Vb/2=QAuo. The solid lines 
represent (7=0, the dashed lines U=QAuj. 

U differs by the driving frequency uj between these two 
figures. The heights of the individual steps in the current 
differ considerably as does the physical reasoning for the 
steps. This is best seen if one looks at the partial currents 
7i, /2, and /a which constitute the Tien-Gordon current. 
For the large electron interaction case the current /2, for 
example, sets in much earlier. The Coulombic staircases 
within the partial currents /„ have of course a width de- 
pending on the size of U . Similar studies can be performed 
for an attractive interaction J7 < 0. 

Laser pulse. — Above we investigated laser fields with 
constant amplitude Eq = 2.405w, i.e. equal to a zero of the 
Bessel function Jq. One can also use other external fields 
with time-dependent envelope function, e.g. with a Gaus- 
sian envelope Eo{t) = 2.405a; exp (-(t - tof/{2a'^)) [20]. 
For the following example we set to =400 fs and cr=80 
fs. Since the amplitude of this laser is time-dependent, 
it induces many more inelastic channel states. During the 
pulse the laser amplitude changes Eq from to 2.405u; fol- 
lowing the Gaussian shape. At each moment in time there 
is one resonance state and many PAT-induced energy 
states eoinw with probability J'^{Eo/uj). For a small bias 
voltage only the partial current /q is non-zero and all other 
partial currents vanish. So in this case the fulfillment of 
the CDT condition leads to a complete suppression of the 
current for t — to. For larger bias, the CDT condition 
just suppresses the partial current Iq but not the total 
current. In Fig. 4, when the bias is small, CDT can be 
seen clearly with and without electron interaction. When 
the bias is big, CDT is only visible as a small dip in the 
current. For the case with electron interaction U = 0.4a;, 
the energy level eo + U just lines up with the Fermi level 
for the bias O.So;, and the current timneling through this 



channel is smaller than the resonance current. This is the 
reason why the current is somewhat smaller than without 
electron interaction. 

As mentioned above, the average current / was obtained 
by averaging I(t) over five periods of the highly oscillating 
carrier field. For the laser pulse with time-dependent en- 
velope the current / is of course depending on the number 
Uav of periods used in the averaging procedure. In these 
cases any averaging procedure will not only average over 
the highly-oscillating carrier field but also smooth the en- 
velope function. This effect will become more pronounced 
for large Uav For a Gaussian shape of the laser pulse the 
averaged current will still be close to a Gaussian form al- 
though slightly shifted and unsymmetric. For the example 
in Fig. 4 there are only minor changes when changing Uav 
between 2 and 10. The averaging process does not depend 
much on Uav for monochromatic laser pulses. 

Conclusions. — To conclude, we have shown that the 
phenomenon of CDT already appears in molecular junc- 
tions which can be modeled by one site coupled to two 
contacts and not only in systems in which the tunneling is 
quenched within the molecule itself. Since the total cur- 
rent consists of several partial currents additional channels 
might be opened by the external field. A complete sup- 
pression of the current occurs only for small bias voltages 
when only the channel for Iq is open and if the amplitude 
of the external field fulfills the CDT condition. In the 
present letter we concentrated on the case of a high carrier 
frequency but as in the case of the two-level system similar 
effects should be possible for smaller frequencies [33] . The 
density matrix formalism used here restricts the present 
results to weak coupling between the leads and the molec- 
ular junction. In the case of CDT between two sites the 
phenomenon does not depend on the size of the coupling. 
Therefore we believe that the effect of CDT between a sin- 
gle site coupled to one or more leads will also survive for 
stronger coupling between site and leads. Although the 
present model calculations certainly can only be a very 
rough description of real systems, the hope is that the ba- 
sic physics survives in more realistic and complex systems. 
Then the effect of CDT with short laser pulses might allow 
for the construction of fast opto-electronic switches if one 
finds materials with long enough coherence times. 

The authors would like to thank the DFG for financial 
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